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Abstract 

Using the Langevin approach and the multiscale technique, a kinetic 
theory of the time and space nonlocal fluctuations in the collisional plasma 
is constructed. In local equilibrium a generalized version of the Callen- 
Welton theorem is derived. It is shown that not only the dissipation 
but also the time and space derivatives of the dispersion determine the 
amplitude and the width of the spectrum lines of the electrostatic field 
fluctuations, as well as the form factor. There appear significant differ- 
ences with respect to the non-uniform plasma. In the kinetic regime the 
form factor is more sensible to space gradient than the spectral function 
of the electrostatic field fluctuations. As a result of the inhomogeneity, 
these proprieties became asymmetric with respect to the inversion of the 
frequency sign. The differences in amplitude of peaks could become a new 
tool to diagnose slow time and space variations in the plasma. 
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Fluctuations find an application in diagnostic procedures. Indeed, 
plasma parameters such as temperature, mean velocity, density and 
their respective profiles can be determined by incoherent (Thom- 
son) scattering diagnostics pQ, i.e. by the proper interpretation of 
data obtained from the scattering of a given electromagnetic field 
interacting with the system. The key point of interpretating them 
is the knowledge of the intensity of the dielectric function fluctu- 
ations or equally of the electron form factor (5n e 5n e )ujs.- Here ui 
and k are respectively the frequency and wavevector of the autocor- 
relations. Due to the Poisson equation the electron form factor in 
the spatially homogeneous system is directly linked to the electro- 
static field fluctuations, which have been the object of active study 
since the early 1960s £Q. In the thermodynamic equilibrium, the 
electrostatic field fluctuations satisfy the famous Callen-Welton 
fluctuation-dissipation theorem [2]: 



(5meu = e 8 ^;? (i) 



1 



linking their intensity to the imaginary part of the dielectric function e{u),\s.), 
and the temperature in energy units. The spectral function Q has peaks, 
corresponding to proper plasma frequencies. The matter becomes more tricky in 
the non-equilibrium case, when the state of the plasma is given by Maxwellian 
distributions characterized by different constant temperatures and velocities 
per species (0 Q , V a ; a = e, i). We have indeed shown [3], that, in the collisional 
regime equations should be revisited. We stressed the fact that a kinetic 
approach should be taken. Introducing fluctuations by the Langevin method, 
we have elaborated a "revisited" Callen-Welton formula containing, beside the 
terms appearing in Eq. new terms explicitly displaying dissipative non 

equilibrium contributions. 

(<SE<SE) uk = V ^[Jmx«+MQa-Q0*i+M k -Va-k-V 6 )*a], 

a ,b=e. t (w-k.V )|e(w,k)| 

(2) 

where Xa (a = e, i) is the complex dielectric susceptibility of the a-th component 
It is important that these new terms contain the interparticle collision frequency 
v a b, the differences in temperatures (O a — Ob) and velocities (V a — V&), and 
the functions $i and <I>2 of real parts of the dielectric susceptibilities. It is 
however not evident that the plasma parameters - temperature, velocities and 
densities can be kept constant . Inhomogeneities in space and time of these 
quantities will certainly also contribute to the fluctuations. Obviously, to treat 
the problem, a kinetic approach is required , especially when the wavelength 
of the fluctuations is larger than the Debye wavelength. To derive nonlocal 
expressions for the spectral function of the electrostatic field fluctuation and 
for the electron form factor we use the Langevin approach to describe kinetic 
fluctuations 0J|S]. The starting point of our procedure is the same as in [3J. A 
kinetic equation for the fluctuation Sf a of the one-particle distribution function 
(DF) with respect to the reference state f a is considered. In the general case the 
reference state is a none-equilibrium DF which varies in space and time both 
on the kinetic scale ( mean free path l e i and interparticle collision time Vei 1 ) 
and on the larger hydrodynamic scales. These scales are much larger than 
the characteristic fluctuation time In the non-equilibrium case we can 

therefore introduce a small parameter fi = v e i/uj, which allows us to describe 
fluctuations on the basis of a multiple space and time scale analysis. Obviously, 
the fluctuations vary on both the "fast" (r, t) and the "slow" {fir, fit) time and 
space scales: 5f a (x,t) — 5f a (x,t, fit, fir) and / a (x,t) = f a {p, fit, fir). Here x 
stands for the phase-space coordinates (r, p). The Langevin kinetic equation 
for df a has the form 0J|3| 

W<Va(x,i) - */£ (x,i)) = -e a SE{r,t)-^^-, (3) 

where 

La* = f + v • £ + f a {t, r); f tt (t, r, p) = e a E-^ - SI a ; 
5I a is the linearized Balescu-Lenard collision operator. 
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The Langevin source in Eq. @ is determined |3] by following equation: 

L axt 5f a (xJ)Sf b (x>,t>f = S ab 5{t - t')S(x - x')/„(x',i'). 
The solution of Eq. ((3J) has the form 

t 

5f a ( X ,t)=5f s (x,t)~J2 J dx ' / dt'Gabfct, x'/ME(r', 1') 91 ^^ , (4) 

k — OO 

where the Green function G a &(x,t, x',t') of the operator L ax ^ is determined by 

£ axt G afc (x,i, x',i') = * o6 <J(x - x')<y(t - «') 

with the causality condition G a f,(x,t,x',i / ) = 0, wheni < t'. Thus, 5f a (x,t)5fb(x.',t') 
and G a h(x,t,x',t') are connected by the relation: 

5/a(x,t)J/ 6 ( X ',i') S = G ah (x,t,x',t')/b(x'/)- 

For the stationary and spatially uniform systems, when DF / a and the op- 
erator T a do not depend on time and space, G a ;,(x,i, x',i') can depend only on 
its time and space variables through the difference t — t' and r — r'. In the 
general case, when the one-particle DF f a (p,fir,fit) and the operator T a slowly 
(in comparison with the correlation scales) vary in time and space, and when 
non-local effects are considered, the time and space dependence of G ab (x,t, x',i') 
is more subtle. 



G a6 (x,i, x'/) = G ab (p, p', r - r',t-t', fir', fit'). (5) 

For the homogeneous case this non-trivial result was obtained for the first 

time in For inhomogeneous systems it has been generalized recently in [7|. 

The relationship JSJ is directly linked with the constitutive relation between 

the electric displacement and the electric field: 
t 

D i (r,t)=jdr' f dt'e ij (r,r',t,t')E j (r',t'). 

— oo 

Previously two kinds of constitutive relations were proposed phenomenolog- 
ically for a weakly-inhomogeneous and slowly time-varying medium: 

(i) the so-called symmetrized constitutive relation 8 : 

A(r,i) = Jdr'J J ( ft eij .(r-r / ,t-t , ;Ai^L,M^y-))Sj(r',*')- (6) 

— oo 

(ii) the non- symmetrized constitutive relation 9 : 

t 

A(r,*) = J dr' J J dtSijir-r'j-t'ifir^tyEjir^t'). (7) 

— oo 

Both phenomenological formulations (i) and (ii) are unsatisfactory. The correct 
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expression should be 

t 

A(r,t) = J dr' J J dte^(r-r',t-f / ;/ir',/ii / ))-E J (r , ,t')- (8) 

— OO 

Taking into account the first-order terms with respect to fi from and JSJ 
we have 

oo 

Sf a (x,t) = <J/£(x,f) d P' d P J dr 



^ 1T ~d(A ~^ P ' £)p~x ' eb<SE ( r_ ' ' t-vGabip, t, p, p , /Ltt, ^r) — , (9) 

(p = r-r',r = t-t'). 

^From the Poisson equation 

^ E(r ' t} = -^E e " / T^Wx', (10) 
and performing the Fourier-Laplace transformation 

oo 

<5E(k, w) = / eft / dr6E(r, f) exp(-At + iwt - ik • r). 

o 

from © we have 

<5E(k,w, fit, fir) = 5W(k,uj) + 4 ™ e a J d p[{l + i & 



dui dfit 



k 7-i jst^a. .. ..^ df a (p,fir,fit) 



-i—SE^.rt-^ "—-' I. (11) 

Here and in the following for simplicity we omit /i, keeping in mind that deriva- 
tives over coordinates and time are taken with respect to the slowly varying 
variables. The resolvent in Ijllfl is determined by the following relation: 

oo ^ 

J dp J dr exp(-Ar + iujr - ik ■ p)G ab (p, r, p, p',fit, fir) = L~^ k d ab S(p - p'). 
o 

The approximation in which Eq. (|llfl was derived corresponds to the geo- 
metric optics approximation 10 . At first-order and after some manipulations, 
one obtains from Eq. ifTTl) the transport equation in the geometric optics ap- 
proximation, which is not considered in the present article, and the equation for 
the spectral function of the electrostatic field fluctuations: 

Ree(io, k) [(6E6E) u , k - - \ (<JE<JE)* k ] =0, (12) 
|e(w,k)| 
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where we introduced 

e(uj,k) = 1 + J2Xa(uJ,k); e(w,k) = 1 + EXa(w,k) 

a a 

^.k) = (l+^|-^-^)^ ) k, i) r) ) (13) 

and where 

Xo(w, k, t, r) = / dpL~j k k • £/ (p,t, r) 

is the susceptibility for a collisional plasma. In the same approximation the 
spectral function of the Langevin source (<5E#E)~, k takes the form 

(5E<5E)5, k = 32 7 r 2 ^e^e J dp (l+i^-i^A)^L-J k f a (p,r,t). (14) 

If Ree(oj, k) ^ 0, it follows from Eqs. ((12(1 and l|14fl that the spectral func- 
tion of the nonequilibrium electrostatic field fluctuations is determined by the 
expression: 



32^ 2 E elRe J dp(l + i& « - i£ • & )£L^ k /„(p, r, i) 
[(<5E*Ekk = 2 — -77 [3 ■ (15) 

The effective dielectric function e(uj, k) in the denominator of Eq. 115(1 de- 
termines the spectral properties of the electrostatic field fluctuations and its 
imaginary part 



_ d d d d 

Ime(u), k) = Jme(u, k) + ——Ree{uj, k, t, r) - — • —Ree(u, k, t, r), (16) 

determines the width of the spectral lines near the resonance. Note that when 
expanding the Green function in Eq. in terms of the small parameter \i, 
there appear additional terms at first order. It is important to note that the 
imaginary part of the dielectric susceptibility is now replaced by the real part, 
which is greater than imaginary part by the factor fi . Therefore, the second 
and third terms in Eq. 1(16(1 in the kinetic regime have an effect comparable to 
that of the first term. At second order in the expansion in /i the corrections 
appear only in the imaginary part of the susceptibility, and they can reasonably 
be neglected. It is therefore sufficient to retain the first order corrections to solve 
the problem. 

For the local equilibrium case where the reference state f a is Maxwellian, we 
have the identity: / dp(l+z£f -z J^^L-J^p, i, r) = JL J dp/ B (p, t, r)- 

— Xa{w, k) (uj a = uj — kV a ) and Eq. 1(15(1 takes the form 

= J2 ** 6 ° 2 Irnxa(u, k). (17) 
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In this case the small parameter /i is determined on the slower hydrodynamic 
scale. For the case of equal temperatures and V a = one obtains a generalized 
expression for the Callen-Welton formula: 

(SESE) UM = 87r6 ^ k) - (18) 
to |e(w,k)| 

To calculate explicitly (JEJE)^ we will restrict our analysis to the vicinity 
of the resonance, i.e. to = ±loq, where Ree(coo, k) = 0. We can develop s(uj, k) = 
(w - UQsgn^Q^lu^sgnu +i[Ime + (JUj - J- • ^)E£e][u=uoBgnu- Thus 



7 8ttT I 

(oj— u>Qsgnu>)' 2 +7 2 ujdRee / du 



(<5E£E) w , k = ,„J„^ , ^J/^ L^o, where 



d 2 d d dRee 

7 = [Ime + dUdt Re£ ~ dk ■ flr ^~flw~ (19) 



is the effective damping decrement. For the case where the system parameters 
are homogeneous in space but vary in time, the correction is still symmetric 
with respect to the change of sign of to, but the intensities and broadening are 
different, and the intensity integrated over the frequencies remains the same as in 
the stationary case. However, when the plasma parameters are space dependent 
this symmetry is lost. The spectral asymmetry is related to the appearance of 
space anisotropy in inhomogeneous systems. The real part of the susceptibility 
Ree is an even function of to. This property implies that the contribution of the 
third term to the expression of the damping decrement (|19[) is an odd function 
of lo. Moreover this term gives rise to an anisotropy in k space. 

Let us estimate this correction for the plasma mode (too = to£) Ree = 
1 _ 4(1 + 3-^), Ime =4^, w? = = 2& and 

~ r 2 dn ljt dn . , , „ 

7 =K i + -- + 6^k.- S ^]/2. (20) 

For the spatially homogeneous case there is no difference between the spectral 
properties of the longitudinal electric field and of the electron density. They are 
connected by the Poisson equation. This statement is no longer valid when 
considering an inhomogeneous plasma. Indeed the longitudinal electric field 
is linked to the particle density by the nonlocal Poisson relation (flop . In the 
latter case, an analysis similar to that made above can also be performed for 
the particle density. ^From Eq. Q there follows 

x n 4\ x s n .,. , ^ 47rike 6 e f . d d , , <9/ a (p,r,i) 
6n a (k,w,r,t)=6nZ(k,w,r,t)+2^ p J dp[(l+t——)L aujk 6n b (k,u,r,t) — 



9 x n 4-\ 9 t-i df a (p,r,t) 

i —Sn b (^,r,t)—L au)k ^ ]. (21) 

At the first order approximation and after some manipulations, one obtains the 
following expression for the electron form factor for a two-component (a = e, i) 
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plasma: 



2n P k 2 



1 + Xi(w,k) 



e(u, k) 



-Jmxe(w,k) 



Xe(w,k) 



E(w,k) 



Qi 2n e k 2 



7mxi(w,k), 



(22) 



where we used for local equilibrium the following expression for the "source" 
(Sn a 5n b )S k = S ai 



ande(w,k) = l+£ ^(w, k); ^(w, k) = (l+i-^^-i^-^-kj-^-k^XaiuJ, k, i, r). 

a 

As above we can expand e(cj, k) near the plasma resonance u) — lol- Thus, for 
the electron line, 



(£rc e £ne) w ,k 
where 



2n.fc^ 



(w-si3nwi,) 2 +(^) 2 uk^dRee/du L^-wz, ) 



7 = [/me + 



<9 2 i?ee 1 9 
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kjRee]/ 



dRee 



uj—ULSgnuj 



(23) 



is the effective damping decrement for the electron form factor. At this stage 
of calculation, let us note that the damping decrements for the electrostatic 
field fluctuations [Eq. (|19() ] and for the electron density fluctuations [Eq. (|23() ] 
are not the same. The origin of this difference is that the Green function for 
electrostatic field fluctuation and density particle fluctuations are not the same. 
This property holds only in the inhomogeneous situation. An estimation for the 
plasma mode is then: 
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2 dn 
n dt 



6k 



2 )sgnu]/2. 



(24) 



^From this equation we see that the inhomogeneous correction in Ea. (|24|l is 
greater than the one in Eq. I|2(JI) by the factor 1 + k 2 D /6k 2 . For the same 
inhomogeneity; i.e., the same gradient of the density, we plot the form fac- 
tor (Sn e 5n e ) u ^ together with the (JEiJE)^ as functions of frequency (Fig. 
1). This figure shows that the asymmetry of the spectral lines is present both 
for (Sn e 6n e )ijj ^ and (^E^E)^. However, this effect is more pronounced in 
(6n e 6n e )u,k than in (<5E(5E) Wj k- 



Conclusion 1 We have shown that the amplitude and the width of the spectral 
lines of the electrostatic field fluctuations and form factor are affected by new 
non-local dispersive terms. They are not related to Joule dissipation and appear 
because of an additional phase shift between the vectors of induction and electric 
field. This phase shift results from the finite time needed to set the polarization 
in the plasma with dispersion. Such a phase shift in the plasma with space dis- 
persion appears due to the medium inhomogeneity. These results are important 
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for the understanding and the classification of the various phenomena that may 
be observed in applications; in particular, the asymmetry of lines can be used as 
a diagnostic tool to measure local gradients in the plasma. 
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Figure 1: The electron form factor ((>n e c>n e ) Wi ic ( solid line) and the spectral 
function of electrostatic field fluctuations (5E5E) w ,k (dashed line) as a function 

of frequency, k-f* = k -f = 6 
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